CALT-68-2110 

JHU-TIPAC-97009 

hep-ph/9705235 



Vub from the Hadronic Invariant Mass Spectrum 
in Semileptonic B Decay 



Adam F. Falk'^, Zoltan Ligeti'' and Mark B. Wise^ 

"'Department of Physics and Astronomy, The Johns Hopkins University 
3400 North Charles Street, Baltimore, MD 21218 

^California Institute of Technology, Pasadena, CA 91125 



Abstract 

The hadronic invariant mass spectrum for the inclusive charmless semileptonic 
decay B e is studied. Particular attention is paid to the region su < 

nn?£), which may be useful for extracting the value of \Vub\- The sensitivity 
of the spectrum to the parameter A = — is explored. Perturbative 
QCD corrections to dV/dsH of order a^/3o are calculated. For sh ~ Am^ 
nonperturbative QCD effects are important and the shape of the invariant 
mass spectrum is controlled by the B meson matrix element of an infinite sum 
of local operators. The utility of the hadronic mass spectrum for extracting 
\Vub\ is explored. 
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The traditional method for extracting \Vub\ from experimental data involves a study 
of the electron energy spectrum in inclusive charmless semileptonic B decay For a 
particular hadronic final state X the maximum electron energy is 

(in the B rest frame), and consequently electrons with energies in the endpoint region 
ms/2 > Ef. > {m\ — m|))/2ms (neglecting the pion mass) must arise from the h ^ u 
transition. A determination of \Vub\ from experimental data on the electron spectrum in the 
endpoint region is possible, provided a theoretical prediction for the electron spectrum can 
be made. 

Recently there has been considerable theoretical progress in our understanding of inclu- 
sive semileptonic B decay [0,^3 . It is based on the use of the operator product expansion 
(OPE) and heavy quark effective theory (HQET) to include in the differential decay rate 
nonperturbative effects suppressed by powers of Aqcd/'^6- At leading order in the Aqcd/''^6 
expansion the B meson decay rate is equal to the h quark decay rate. There are no non- 
perturbative corrections of order AqcD/^^b- At order Aqcp/m^ [|],|[ the nonperturbative 
corrections are characterized by two HQET matrix elements Ai^2; which are defined by 

Ai = (i?(t;)|/i(^) {iDfh^'^\B{v))/2mB, 

A2 = {B{v)\ I h^^ a,,G^^ h^^ \B{v))/QmB . (1) 
These matrix elements also occur in the expansion of the B and B* masses in powers of 

niB = mf, + A - (Ai + 3A2)/2mfe + . . . , 

niB* = mf, + A - (Ai - A2)/2mf, + . . . . (2) 

Similar formulae hold for the D and D* masses. The parameters Ai and A2 are independent 
of the heavy h quark mass (there is a weak logarithmic dependence in A2) and are of order 
Aqcd- The measured B* — B mass splitting fixes A2 = 0.12 GeVl The mass formulae for 
the B and B* mesons involve not only Ai,2 but also a parameter A, which is the difference 
between the B meson mass and the b quark mass in the —>■ 00 limit. The measured 
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B semileptonic decay spectrum in the region > 1.5 GeV has been used to determine 
A ~ 0.4 GeV and Ai — 0.2GeV^ [^]. Unfortunately the uncertainties from terms of order 
{-^QCB / f^b)^ are quite large 0. (A linear combination of A and Ai is rather well constrained, 
but the individual values are more uncertain.) 

The maximum electron energy in semileptonic b quark decay is mf,/2. This is less than the 
physical endpoint by A/2, which is comparable in size to the endpoint region Ai^^^'^'^P"™*-' = 
m|)/2mB ^ 0.33 GeV. Using the operator product expansion and HQET, the effects which 
extend the electron spectrum beyond its partonic value appear as singular terms in the 
prediction for dT/dEe involving derivatives of delta functions, S^'^\Ee — mb/2). Near the 
endpoint the electron spectrum must be smeared over a region of energies AE^ before theory 
can be compared with experiment. If the smearing region AEf, is much smaller than Aqcd, 
then higher dimension operators in the OPE become successively more important and the 
OPE is not useful for describing the electron energy spectrum. For AEg much greater 
than Aqcd; higher dimension operators become successively less important and a useful 
prediction for the electron spectrum can be made using the first few terms in the OPE. 
When AEe ~ Aqcd there is an infinite series of terms in the OPE which are all equally 
important. Since AE^^^^p™'^*) is about Aqcd, it seems unlikely that predictions based on a 
few low dimension operators in the OPE can successfully determine the electron spectrum 
in this region. 

In the future, another possibility for determining \Vub\ may come from a comparison of 
the measured hadronic invariant mass spectrum in the region sh < fn'jj with theoretical pre- 
dictions. Here Sh = {pb — qY, where is the B meson four-momentum, and q = Pe+Pue is 
the sum of the lepton four-momenta. An obvious advantage to studying this quantity rather 
than the lepton energy spectrum is that most of the B —>■ e Ue decays are expected to lie 
in the region sh < rnj^, while only a small fraction of the B —>■ e Vf. decays have electron 
energies in the endpoint region. Both the invariant mass region, sh < rnj^, and the elec- 
tron endpoint region, 771^/2 > E^ > {Tn\ — mjj)/2mB, receive contributions from hadronic 
final states with invariant masses that range up to m^i. However, for the electron endpoint 
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region the contribution of the states with masses nearer to rriD is kinematically suppressed 
since they typically decay to lower energy electrons. In fact, in the ISGW model [j^ the 
electron endpoint region is dominated by the vr and the p, with higher mass states making 
only a small contribution. The situation is very different for the low invariant mass region, 
sh < ^"d, with no cut on the electron energy. Now all states with invariant masses up 
to niD contribute without any preferential weighting towards the lowest mass ones. In the 
ISGW model the vr and the p mesons comprise only about a quarter of the B semileptonic 
decays to states with masses less than rriD. Consequently, it is much more likely that the 
first few terms in the OPE will provide an accurate description B semileptonic decay in the 
region Sh < tn'jj than in the endpoint region of the electron energy spectrum. Combining 
the invariant mass constraint, sh < ^"d, with a modest cut on the electron energy will not 
destroy this conclusion. (Such a cut will probably be required experimentally for the direct 
measurement of sh via the neutrino reconstruction technique.) We also expect that the 
B Xu e z/g rate in the invariant mass region sh < rnj^ is less sensitive to nonperturbative 
effects than is the rate in the hadron energy region Eh < ttlj:, (in the B rest frame) [Q, 
since the hadron energy constraint cuts out more of the phase space for states with mass 
near m/5 than for the lower mass states. In this letter we explore the utility of the hadronic 
invariant mass spectrum for determining the magnitude of Vub- The possibility of using 
the hadronic invariant mass spectrum in B —>■ XcCUe to determine A and Ai was discussed 



in Ref. |T^. The technique is promising but awaits better data on dT/dsn- 

To begin with, consider the contribution of dimension three operators in the OPE to the 
hadronic mass squared spectrum in B X„ e Ue decay. This is equivalent to b quark decay 
and implies a result for dT/dEodso (where Eq = pf, ■ {p^ — q)/mb and Sq = {pb — qY are 
the energy and invariant mass of the strongly interacting partons arising from the h quark 
decay) that can easily be calculated using perturbative QCD up to order a1(3o. Even at this 
leading order in the OPE there are important nonperturbative effects that come from the 
relation between the h quark mass and the B meson mass in Eqs. (0). The most significant 
effect comes from A, and including only it {i.e., neglecting the effect of Ai^2); the hadronic 



invariant mass sh is related to sq and Eq by |T0 



sh = so + 2AEo + . (3) 

Changing variables from (sq, Eq) to {sh, Eq) and integrating Eq over the range 

Vs^-A<^o < 77^(sH-2AmB + m|), (4) 

gives dr/dsH, where A^ < < ""^b- Feynman diagrams with only a w-quark in the final 
state contribute at sq = 0, which corresponds to the region A"^ < sh < Attib- 

Although dT/dsH is integrable in perturbation theory, it has a double logarithmic sin- 
gularity aX Sh = Atub- At higher orders in perturbation theory, increasing powers of 
as\v?[[sH — AniB) /rn?B\ appear in the invariant mass spectrum.B Therefore, dV/dsu in 
the vicinity of sh = Atub is hard to predict reliably even in perturbation theory. (In the 
region sh ^ Atub nonperturbative effects, which we discuss later, are also important.) The 
behavior of the spectrum near sh = Auib becomes less important for observables that av- 
erage over larger regions of the spectrum, such as dV/dsu integrated over sh < A^, with 

significantly greater than Aitlb- Therefore, we present results for dF/ds/^ in the region 
Sh > AniB, where only the bremsstrahlung Feynman diagrams contribute. Calculating these 
Feynman diagrams gives the differential decay rate 

^2 ™3 / A \ 3 



dsH 1927r3 V tub) 



X 



X{sH,A) + i f3QY{sH,A) + ... 



2 



(5) 



where /3o = 11 — 2r2j/3 is the one-loop beta function of QCD. 

In Figs. 1 we plot X{sh,A) and Y{sh, A) as functions of sh for A = 0.2, 0.4 and 0.6 GeV. 
The MS scheme is used for the strong coupling, and we choose to evaluate as at the scale 
^/sH. While Y{sh,A) is sensitive to this choice, the sum of the two terms in the square 
brackets in Eq. has only a weak scale-dependence. Even though the a^jSQ correction is 



*For recent discussions of a similar phenomenon in the electron energy spectrum, see Ref. |11| 
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FIG. 1. The functions X{sh,A) and Y{sh,A) defined in Eq. (|) for A = 0.2 GeV (dotted 
curve), 0.4 GeV (solid curve), and 0.6 GeV (dashed curve). 

as large as the «<, term, this does not necessarily imply a problem with the perturbative 
corrections, since there is a renormalon ambiguity of order Aqcd in A which cancels a 
renormalon ambiguity in the perturbative QCD corrections. 

To examine the sensitivity to A of an extracted value of \Vub\ from the number of events 
in a region s// < A^, we define the dimensionless quantity f (A^, A) by 

.5 / X . 5_ 



dsH 1927r3 



1-— ) r(A^A). 

UIbJ 



(6) 



In Fig. 2 we plot f (A^,A) as a function of A^ for A = 0.2, 0.4 and 0.6 GeV in the region 
KniB < A^ < 4.5 GeV^, using ^^(mb) = 0.2. These curves approach r(m^,A) ~ 0.73 



as A^ 



m 



B 



[121. The spread of the curves in Fig. 2 together with the (1 — A/itlb)^ 



dependence factored out in Eq. (||) suggest that an accurate value of \Vub\ can be obtained 
from the number of events in a region sh < A^ if A^ is not much below m|, = 3.5 GeV^, and 
if a reasonably precise determination of A is available. For example, with A^ = 3.5 GeV^ 
and A = 0.4 ±0.1 GeV, the uncertainty arising from the error in A in the extracted value of 
\Vub\ is only 8%. So far nonperturbative corrections from higher dimension operators in the 
OPE have been neglected. We discuss their influence on the extraction of \Vub\ later. 
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FIG. 2. The function f (A^, A) defined in Eq. (|) as a function of for A = 0.2 GeV (dotted 
curve), 0.4 GeV (sofid curve), and 0.6 GeV (dashed curve). 

Experimental uncertainties will cause some of the B ^ XcCVe events to appear to have 
sh < fnj^. If experimental sh resolution forces to be much below m|,, the uncertainties 
increase significantly. It is not clear at the present time how, for example, = (1.5 GeV)^ 
compares to Am^. For such a small value of A^, our results for f(A^,A) in Fig. 2 are 
only reliable if A has a small value, below 0.4 GeV. With A = 0.4 GeV, one should worry 
about the reliability of an extraction of \Vub\ based on A^ = (1.5 GeV)^, since higher order 
perturbative corrections and nonperturbative effects (which we discuss next) are likely to 
be important.! 

In the low mass region Sh ~ ^itib, nonperturbative corrections from higher dimension 



'''For example, the order a'^Po result predicts for A = 0.4 GeV that a large fraction (about 40%) 
of the B XuCUe events have sh > (1.5 GeV)^. Taking Fig. 2 literally and assuming A = 
0.3 lb 0.1 GeV, the uncertainty in \Vub\ would be 17%, but the sensitivity to uncalculated higher 
order perturbative and nonperturbative effects could be significant. 
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operators in the OPE are very important. Just as in the case of the electron spectrum in 
the endpoint region , the most singular terms can be identified and summed into a shape 
function, S{sh)- Neglecting perturbative QCD corrections, we write 



dsH 1927r3 

It is convenient to introduce the scaled variable y = sn/^iTih and define a dimensionless 
shape function S{y) = AmbS{sH)- Then 

siy) = T^y^r^n ^ k"^^ (3 - 2y) e{i - y)] . (8) 

The matrix elements An are the same ones that determine the shape functions for the 
semileptonic B decay electron energy spectrum in the endpoint region and the endpoint 
photon energy region in weak radiative B decay. Explicitly, 

{B{v) I h^^ iD^, . . . iD,^ h^^^ \B{v))/2mB = AnV,,... v,^ 

+ terms involving the metric tensor . (9) 

The AnS have dimension of [mass]", and hence the coefficients A„/A" are dimensionless 
numbers of order one. The first few A^s are Ao = 1, Ai = 0, A2 = — Ai/3, and ^3 = — pi/3. 
Using the equations of motion, pi can be related to the matrix element of a four-quark oper- 
ator. In the vacuum saturation approximation, pi = (27r 0:5/9)7715/^ Unfortunately, 
the scale-dependence of this result leaves the value of pi highly uncertain 0. 

The shape function S{y) is an infinite sum of singular terms which gives an invariant 
mass spectrum that leaks out beyond y = I {i.e., sh = Ami,). For y ~ 1 {i.e., sh ~ Am;,) all 
terms in Eq. (^) are formally of equal importance. Since Am^ ^ 2 GeV^ is not too far from 
m^, it is necessary to estimate the infiuence of the nonperturbative effects on the fraction 
of B decays with invariant hadronic mass squared less than m|). It is difficult to obtain a 
model-independent estimate of the leakage of events above an experimental cutoff s// = A^, 
given that we can estimate only the first few moments. An. Upper bounds on this leakage 
can be obtained if S{y) is assumed to be positive definite. This is consistent with the naive 
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interpretation of the leading singularities as constituting a nonperturbative smearing of the 
rate beyond the b quark decay endpoint |T^; however, there is no proof that this property 
actually holds. The differential decay rate is positive, but for sh comparable with there 
are other nonperturbative terms which are equally important. Furthermore, perturbative 
QCD corrections have been neglected. Models, such as the ACCMM model |TB[, do give a 
positive shape function. 

The fraction of events with sh < is given by 

F{A)= r^'^^ySiy), (10) 
Jo 

where e(A) = A^/Am;,. Recall that the kinematic point sh = A^ corresponds to y = e. 
Assuming a positive shape function, F{A) is greater than 

Fp(A)=/ dyP{y,eiA))Siy), (11) 
Jo 

provided P{y,e) satisfies the following properties: (i) P{y,^) < 1 for ?/ < e; (m) P{y,e) < 
for y > e. The lower bound F{A) > Fp{A) holds for any such P{y,e). Furthermore, if 
P{y,e) = Pk{y,e) is a polynomial of degree k in y, then only the first k moments, Ak, 
appear in the bound. Setting Pkiy, e) = Z^o (^ei^) y^i ^"^^ integrating by parts n times yields 

As an illustration of the utility of this bound, consider first the simple quadratic poly- 
nomial P2{y) = 1 — y^ 1^ ■ This leads to the lower bound 

F(A)>l-^ + ^, e = AVAm,. (13) 

For A = 0.4 GeV and Ai = -0.2 GeV^ the bound is Fimo) > 76%. With larger A, the 
bound weakens dramatically. For A = 0.6 GeV and Ai = — 0.2GeV^, it is only F^mc) > 
59%. Once again, an independent determination of A and Ai is necessary for these bounds 
to become useful. For a cubic polynomial, we also need to know pi. For example, consider 
P^iy) = 1 ~ y^ 1^ ■ Then the lower bound is 
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^(^)>l-|3+^ + ^' e = AVAm.. (14) 

For A = 0.4 GeV and Ai = -0.2 GeV^, F(m^) > 83% if pi = and F(mz)) > 87% if 
Pi = 0.1 GeV'^. The dependence on A is still the most important, as the latter bound falls 
to F^nio) > 69% for A = 0.6 GeV. We could improve these bounds by optimizing the 
coefficients in the polynomial Pk- Since the optimization itself will depend on A, Ai and pi, 
it does not seem worth while to proceed along this line at the present time. 

If due to experimental resolution one can only use sh < (1.5 GeV)^, then the bounds 
become much weaker. For example, using the cubic polynomial above with A = 0.4 GeV 
and Ai = -0.2 GeV^ the bound is F(1.5 GeV) > 36% for pi = and F(1.5 GeV) > 51% for 
pi = 0.1 GeV^ 

An alternative is to resort to models for an estimate of the effect of high order terms 
in the sum in Eq. (|^). As an example, consider the ACCMM model [|l^,|l^], where the B 



meson is modeled by a spectator quark with mass m^p and momentum p, and a b quark with 
momentum —p and effective mass ml^^^ = tub — J^lp + p"^- The probability that the spec- 



tator quark momentum takes the value p is $(p ). In this model A = / d^p $(p ) \j^1p + P'^^ 
and Ai = — / d^p$(p) \p\^. 

To plot the shape function S{sh) in the ACCMM model, we neglect the boost from the 
h quarkfe quark rest-frame into the B meson rest-frame (such affects are subleading in the 
rrib —>■ oo limit for all values of sh)- In Fig. 3 we plot the shape function for three different 
cases that give A = 0.4 GeV. They are $(p) oc e~'^'/^^ with pp = 0.13 GeV and m^p = 
(dotted curve); $(p) oc e^'^'^/^f' with pp = 0.35 GeV and m^p = (dashed curve); and 
Pf = 0.3 GeV and m^p = 0.2 GeV (solid curve). In these cases only 4.4%, 4.1% and 2.8%, 
respectively, of the B X^eUe decays have sh > m-'j^. Even if the leakage into the region 
Sh > fTiD were a factor of two or three greater than this (a possibility which is not at all 
unlikely), unknown nonperturbative effects characterized by the A„'s only give rise to about 
a 10% uncertainty in the fraction of B semileptonic decays with sh < iTij-,. On the other 
hand, if due to charm contamination one can only use events with sh < (1.5 GeV)^, then 
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FIG. 3. The shape function S{sh) in the ACCMM model, ^{p) cx e'lPl/P^ with pp = 0.13 GeV 
and rrisp = (dotted curve); ^{p) oc e~^^^^^^F with pp = 0.35 GeV and nisp = (dashed curve); 
and pf = 0.3 GeV and m^p = 0.2 GeV (solid curve). 

the sensitivity to the shape function is much greater. For the three models in Fig. 3, the 
fraction of B XuCUe decays with sh > (1.5 GeV)^ is 15%, 16%, and 15%, respectively. 
One should not conclude from the approximate agreement between these models that the 
uncertainty in these predictions for the leakage is less than a factor of two. 

The larger the value of A, the larger the fraction of B ^ X„ e Ue decays that leak out 
beyond sh = fn'j^. In Fig. 4 we plot the model $ oc e^'^'^/^?" with m^p = 0.2 GeV and pp 
taking three values {pp = 0, 0.3, and 0.5 GeV) corresponding to the choices A = 0.2, 0.4 and 
0.6 GeV. For these values of pp, the model gives Ai = 0, —0.14, and —0.38 GeV^ respectively, 
in qualitative agreement with the correlation between A and Ai from the electron energy 
spectrum in semileptonic B decay The pp = (dotted) curve is given analytically 

by the n = term in Eq. (§) with A = 0.2 GeV. The fraction of B e decays with 

Sh > "n^z) is 0, 2.8%, and 12%, respectively. The fraction of events with sh > (1-5 GeV)^ is 
0, 15%, and 31%, respectively. The rapid variation of these values with A shows again that 
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FIG. 4. The shape function S{sh) in the ACCMM model, with m^p = 0.2 GeV. The dotted 

curve corresponds to A = 0.2 GeV (pp = 0), the sohd curve is A = 0.4 GeV {pp = 0.3 GeV), and 
the dashed curve is A = 0.6 GeV {pp = 0.5 GeV). 

a reliable determination of \Vub\ from the number of events in the region sh < is only 
possible if A does not have too large a value. This is especially true if experimental issues 
force to be significantly smaller then m|,. 

We have investigated the utility of the hadronic invariant mass spectrum in B ^ X„ e z/e 
decay in the region sh < rnjj for a possible model-independent determination of \Vub\- 
Perturbative QCD corrections to dV/dsH of order a'^Po were calculated paying particular 
attention to kinematic effects arising from A. A measurement of dsn (dF/dsj/) can be 
translated into a value of \Vub\ using Fig. 2, and the result can then be corrected (with some 
model dependence) for nonperturbative effects coming from operators with dimension five 
and higher in the OPE. If A^ ~ m'j^ and A = 0.4 ±0.1 GeV are experimentally feasible, 
then \Vub\ can be extracted model- independently with about 10% theoretical uncertainty 
from the error in the value of A. In this case uncertainties associated with higher dimension 
operators in the OPE are likely to be small. A determination of A with a precision of 
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±0.1 GeV from experimental information on semileptonic B decay and weak radiative B 
decay seems possible P,P, p!0|JT7| , p!8[] . 

It is possible that the experimental invariant mass resolution will necessitate an upper 
cut on the hadronic invariant mass squared which is somewhat below mj^, in a region 
where nonperturbative effects that make the spectrum leak beyond are not negligible. For 
Sh ~ ^fns, both nonperturbative strong interaction effects and higher order perturbative 
corrections become important. If A^ has to be substantially smaller than mj^, then A 
cannot be too large for this method of extracting \Vub\ to remain viable. In this case, the 
theoretical uncertainty in \ Vub\ will depend sensitively on both the value of A^ and the value 
and uncertainty in A at the time when a measurement of dsn {dT/dsn) is available. In 
addition, since the nonperturbative effects introduce a certain level of mo del- dependence, it 
will be important to compare the extracted value of \Vub\ from the hadronic invariant mass 
spectrum with its value from other determinations, such as from exclusive decays [|19 . 
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